Exact Results for the Crossover from Gaussian to Non-Gaussian Order Parameter 
Fluctuations in Quasi One-Dimensional Electronic Systems 
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The physics of quasi one-dimensional Peierls systems is dominated by order parameter fluctua- 
tions. We present an algorithm which allows for the flrst time to exactly calculate physical properties 
of the electrons gas coupled to classical order parameter fluctuations. The whole range from the 
Gaussian regime dominated by amplitude fluctuations to the non-Gaussian regime dominated by 
phase fluctuations is accessible. Our results provide insight into the 'pseudogap' phenomenon occur- 
ring in underdoped high-Tc superconductors, quasi one-dimensional organic conductors and liquid 
metals. 
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Pseudogaps have been discussed recently in the con- 
text of the high temperature superconductors. Methods 
which were originally designed for studying the quasi 
one-dimensional problems ||l|, ||] have been used to de- 
scribe the physics of the suppression of electronic den- 
sity of states due to antiferromagnetic or superconduct- 
ing fluctuations H, ^, H, ||. Here we present a method 
which allows to describe the crossover from the Gaussian 
regime dominated by amplitude fluctuations to the non- 
Gaussian regime dominated by phase fluctuations with- 
out any approximation. We solve the problem originally 
posed by Lee, Rice and Anderson |l| in 1973 exactly. 

Although the pseudogap problem is more general, it 
first appeared in the context of the charge density wave 
(CDW) systems. The Peierls transition of quasi one- 
dimensional electronic systems like K0.3M0O3 is due to 
the coupling of a particular phonon mode to the elec- 
trons. This coupling would lead to a mean field phase 
transition at some temperature T*^^. However, the 
phase transition of the one-dimensional system to the 
charge density wave phase, which would break a contin- 
uous symmetry is prevented by order parameter fluctua- 
tions. Only at some lower temperature, T^^ , determined 
by the three-dimensional coupling of the one-dimensional 
systems the charge density wave phase is established. 
The Kohn anomaly leads to a softening of the phonon so 
that at some temperature close to the mean field Peierls 
transition temperature, T^^^ the phonon mode can be 
viewed as a static lattice distortion. The properties 
of the electrons in the intermediate temperature regime 
T^'^ > T » T'^^ are determined by the coupling to ID 
CDW order parameter fluctuations, which the mean-field 
theory does not describe even qualitatively. All attempts 
to describe the electronic properties in this regime start- 
ing with Lee, Rice and Anderson ||^, |^, |^, ^, ^ assume 
Gaussian order parameter amplitude fluctuations. Re- 
cent calculations |jlfl, O, 03] corrected a technical mis- 
take in the original paper by Sadovskii [Ej but confirmed 
that the density of states, N{e), of the electrons behaves 
like N{€) ~ e^ below the mean field gap for large cor- 



relation lengths in the Gaussian model which is com- 
pletely unphysical. On the other hand, even a modest 
suppression of the density of states, i.e. a pseudogap, 
requires an enormous correlation length in a Gaussian 
model ||ll|, [1^. Models taking into account phase fluctu- 
ations only ]13|, which should contain the right physics 



far below Tj^ , tend to overestimate the suppression of 
the electronic density at the Fermi surface and cannot de- 
scribe the physics above the mean field transition. Thus 
a more sophisticated approach is needed. 

We begin by defining the problem. The dispersion of 
the electrons close to the Fermi energy can be assumed 
to be linear. The Hamiltonian of the electrons has the 
form: 

H = -ivF{R^d^^R-L^d^L)+A{x)R^L+A*{x)L''R, (1) 

where the operators R^ and L^ create left and right mov- 
ing electrons respectively. The classical order parameter 
field A (a;) is determined by a Ginzburg Landau action 
given below, vp is the Fermi velocity. Contrary to the 
assumption in previous work, it is not sufficient to de- 
scribe the order parameter fluctuations by the variance 
and correlation length only, but that one needs to con- 
sider higher moments of the order parameter correlator. 
This becomes intuitively clear if one considers two cases. 
If the order parameter varies smoothly, as in the Gaus- 
sian regime, regions where the order parameter is sup- 
pressed are smeared out over the correlation length. The 
electronic wavefunction is spread out over a length com- 
parable with the correlation length. The kinetic energy 
is low and consequently many states can be found at low 
energy even when the correlation length is large. On the 
other hand, if the order parameter is established and only 
suppressed over a length scale much shorter than the cor- 
relation length of the potential, as in the non-Gaussian 
regime, the electronic wavefunction decays over a dis- 
tance vp/A and has a large kinetic energy. For the same 
correlation length and variance of the order parameter 
the electronic wavefunction is much stronger suppressed 
for non-Gaussian fluctuations. 



Next we consider the order parameter fluctuations. 
For commensurate fluctuations, the low energy electronic 
density of states is dominated by the Dyson singularity 
which only exists in one dimension |0, n2| . For the more 
general case, the order parameter fluctuations are com- 
plex and the Dyson singularity is absent. Therefore we 
will restrict our discussion to complex order parameters. 
The classical complex order parameter fluctuations, A{x) 
are described by the Ginzburg-Landau functional: 
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(2) 



Close to the mean field phase transition a varies linearly 
with the temperature a{T) — a'{T/T^'^^ — 1), whereas 
b and c (and therefore the length scale, ^o = \/c/a') 
are nearly temperature independent. In principle, the 
coeflficients a, b and c have to be determined self con- 
sistently from the electronic properties. The ID system 
is disordered above the 3D ordering temperature, Tj?^. 
Nevertheless the action Eq. (0) has two different regimes: 
if a{T) is positive and large, the order parameter fluctu- 
ations are centered around zero and basically Gaussian. 
For a{T) negative and large the amplitude of the order 
parameter is given by \/< A^ > and only the phase fluc- 
tuations play a role. 

Here we would like to sum over all configurations 
of the order parameter with the Boltzmann weight, 
exp(— AF[A(x)]//i;Brc)- The technical problem is how 
to generate a sufficiently large configuration of the order 
parameter (typically lengths of a chain: L ~ lO^^o) in 
the intermediate regime so that the electronic properties 
can be calculated reliably. A Monte Carlo simulation of 
such a large system can in principle be done (a more so- 
phisticated algorithm like the Wolff algorithm |1J] has to 
be adopted to avoid critical slowing down close to T^'^) 
but it turns out that there is a much simpler way to per- 
form the calculation. The method presented here is based 
on the transfer matrix formalism first used by Scalapino, 
Sears and Ferrell |jia to calculate the thermodynamic 
properties of classical order parameter fluctuations in one 
dimension exactly. It is useful to write the free energy in 
units where the length is measured in units of ^Qj the 
size of the order parameter in units of Aq = ■y/a'/2& and 



the temperature in units of T^ 
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reduced Ginzburg temperature, Ar, at which the fluctua- 
tions start to dominate is Ar = (a'^/fofesT*^^)"^/'^ (note 
the factor of 2 between our definition and Ref. ||l^ ) . The 
relevant physical parameters are the bare length scale, 
fo, the gap scale, Aq, the mean field critical tempera- 
ture, T^^^, and the size of the fluctuation regime. At. 

The transfer Hamiltonian for Eq. (y) is given in appro- 
priate units by: 
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where A is a two dimensional vector of the real order pa- 
rameter components, A = A' -|- iA" and A = (A', A"). 
The Nabla operator is defined as V = (Sa' , c?A" ) • The 
prefactor of the Hamiltonian, VAr, can be absorbed in 
the length scale. This Schrodinger equation in imaginary 
time is equivalent to a stochastic random walk pq | where 
the groundstate wavefunction is the distribution func- 
tion of the spatial coordinate of the Schrodinger equa- 
tion. The "spatial" coordinate of the transfer Hamilto- 
nian is the order parameter fluctuation A and the imag- 
inary time corresponds to the spatial coordinate along 
the chain x. To see this, we introduce the wavefunction, 
-0, as a ratio '0(A, x) = $(A, a;)/?Ao(A,x) of an auxilary 
function $ and V'o the ground state wave function of the 
anharmonic oscillator for given parameters a, b and c. 
The function $ obeys the following equation of motion: 



— = — V^$ -I- V • — ^ 
dx 2 \ i^o 
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This is nothing but a diffusion equation for $. The dif- 
fusion of the order parameter. A, can thus be described 
by a Langevin equation: 



dA 
dx 
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where fj = {ri',ri"), is uncorrelated Gaussian noise in the 
complex plane which can be generated easily. 

The Langevin equation Eq. (0) is nonlinear but can 
be simulated very easily (see e.g. [0)- The distribution 
of the order parameter fluctuations A is given by con- 
struction by the ground state wave function of the trans- 
fer Hamiltonian. The key observation now is that the 
order parameter fluctuations which can be locally gener- 
ated from the Langevin equation Eq. (m have precisely 
the statistics given by the action Eq. (0). The Langevin 
equation, Eq. (|5|) can be viewed as an extremely efficient 
way to generate typical order parameter fluctuations. 

The method presented here is related to the path in- 
tegral Monte Carlo algorithm in which the solution of 
the Schrodinger equation is obtained by simulating the 
kinetic energy with a diffusion equation and the poten- 
tial energy using a von Neumann rejection to implement 
importance sampling |]l6| . Because the "guiding func- 
tion", ipQ, is the solution of the transfer Hamiltonian 
Schrodinger equation "paths" (configurations) in the or- 
der parameter space are generated according to their 
weight in the partition function. In some sense we have 
inverted the path integral Monte Carlo method to gener- 
ate the paths according to their statistical weight. 

As a first test, we use Eq. (||) to generate a configu- 
ration of a chain with length L = 10^ and calculate the 
average of the square of the order parameter fluctuation, 
< |Ap >, and the correlation length ^ and compare to 
the exact results by Scalapino et al. [|l5| in Fig. |l|. It is 
obvious that even for a relative short chain the calculated 



< 

<f 

A 
V 




LU 




FIG. 1: Variance and correlation length of the order parame- 
ter fluctuations as a function of reduced temperature, r gen- 
erated using Eq. (H) (marked by □) for a finite system with 
length L = 10^ compared to to the exact results by Scalapino 
et al. lis] (dashed line). 



variance and correlation length are basically identical to 
the exact results. The differences are due to the numeri- 
cal evaluation of the drift function V'0o/V-'o- We solve the 
Schrodinger equation Eq. ([l|) numerically and approxi- 
mate the logarithmic derivative of the ground state wave 
function by a fourth order polynom in |A| - which ap- 
parently is a very good representation of |V'0o|/V'o(|^|)- 
As an aside we remark that for the case of Gaussian or- 
der parameter fluctuations the drift term |V?/'o|/''/'o(|A|) 
is linear and corresponds to the Ornstein-Uhlenbeck pro- 
cess considered in Ref. fl^ . 

The density of states of Eq. (|l|) for a given order pa- 
rameter configuration can be calculated by various meth- 
ods, for example by using a lattice version of the Hamil- 
tonian Eq. Im) and exact diagonalization |ll| or some 
more sophisticated method based on the phase formal- 
ism which has been developed recently in Ref. p2| (for a 
detailed description see [Q). We calculate the density of 
states with the Langevin equation for the order parame- 
ter, Eq. (P) which contains amplitude and phase fluctua- 
tions. In this way only local information is needed to cal- 
culate the properties of the electrons. The requirement 
for storage is minimal compared to a full simulation of 
the classical field. The density of states can be obtained 
by differentiating numerically the integrated density of 
states. The resulting density of states as a function of 
energy is shown in Fig. [2[ The en ergy scale at a fixed 
= V< |AP >. 



temperature is A(T) 

The calculated density of states shown in Fig. has 
several features which apparently resolve some of the 
problems of the approximations made in previous cal- 
culations. One of the problems of the Gaussian approxi- 
mations is that only for very large correlation lengths of 
the order parameter fluctuations a pseudogap appears. 
Here we see that already for relatively modest correlation 



FIG. 2: The density of states normalized to the bare den- 
sity of states, A''(0), as a function of energy, E/A{T) and re- 
duced temperature r. The temperature varies from r/Ar = 
— 10 ... in steps of one. As the temperature is decreased be- 
low the mean field transition the density of states is strongly 
suppressed. 



lengths much above the mean field transition the density 
of states at zero frequency is strongly suppressed. Our 
calculation smoothly interpolates between the amplitude 
fluctuation dominated regime and the phase fluctuation 
dominated regime. 
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FIG. 3: The inverse localization length as a function of fre- 
quency for reduced temperature, r varies from r/Ar = 
— 10...0 in steps of one. The localization length increases 
with decreasing temperature. The lower curve corresponds 
to the highest temperature and the topmost to the lowest 
temperature. 

The localization length can be calculated from the 
Thouless relation fia], and is presented in Fig. 0. For 
increasing temperatures the localization length decreases 
and approaches zero uniformly. The localization length 
at a given energy is a monotonic function of tempera- 
ture which is different from the phase fluctuation only 
model [p3. This might be due to the fact that in those 



models the variation of the gap scale as a function of 
temperature is not taken into account. 

Finally we discuss the implication of our calculation for 
experiments. As an example, we calculate the tempera- 
ture dependence of the electronic spin susceptibility as a 
function of temperature. The relevant parameters are the 
size of the gap fluctuations Aq and the size of the fluctua- 
tion regime, At. Apparently excellent agreement can be 
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FIG. 4: Experimentally determined static spin susceptibilities 
Xo{T) for a number of materials from Ref. uM and a tentative 
fit (solid line) using the calculated density of states. The 
parameters are Ao — 2T^ and Ar = 0.1. 



achieved with the experimentally determined Ginzburg 
temperature 20K for K0.3M0O3 |M . At a lower temper- 
ature, TsDj three-dimensional ordering sets in {T^d/Tmf 
= 0.4, 0.6 and 0.26 for TaSg, K0.3M0O3 and (TaSe4)2l re- 
spectively). The departure from the purely ID behavior 
can be most clearly noted for K0.3M0O3. 

To summarize, we have developed a new algorithm 
for studying the electronic properties of the quasi-one- 
dimensional Peierls systems which allowed for the first 
time to study the electronic properties in crossover 
regime from the Gaussian to the non-Gaussian regime of 
the order parameter fluctuations without further approx- 
imations. Similar problems appear in the study of the 
pseudogap phenomenon in the High Temperature Super- 
conductors. More general the question of what happens 
to an electronic system coupled to soft classical degrees 
of freedom is relevant not only for the pseudogap but also 
for the question of how to describe the moment forma- 
tion in an itinerant electronic system. These questions 
have been discussed in terms of sophisticated perturba- 
tion theories, like the Parquet approximation [pl| , or self 
consistent approaches like the SCR theory by Moriya [g2| . 
Our calculation basically demonstrates that it is not pos- 
sible to describe the moment formation by Gaussian fluc- 



tuations (or perturbation theory), but that a self consis- 
tent theory which takes into account the non-Gaussian 
fluctuations into account is necessary. 
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